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ABSTRACT 

We present high-resolution spectroscopy of the neon i^-shell and iron L-shell interstellar absorption 
edges in nine X-ray binaries using the High Energy Transmission Grating Spectrometer (HETGS) 
onboard the Chandra X-ray Observatory. We found that the iron absorption is well fit by an experi- 
mental determination of the cross-section for metallic iron, although with a slight wavelength shift of 
«20 mA. The neon edge region is best fit by a model that includes the neutral neon edge and three 
Gaussian absorption lines. We identify these lines as due to the ls-2p transitions from Ne II, Ne III, 
and Ne IX. As we found in our oxygen edge study, the theoretical predictions for neutral and low- 
ionization lines all require shifts of «20 mA to match our data. Combined with our earlier oxygen edge 
study, we find that a best fit O/Ne ratio of 5.4±1.6, consistent with standard interstellar abundances. 
Our best fit Fe/Ne ratio of 0.20±0.03 is significantly lower than the interstellar value. We attribute 
this difference to iron depletion into dust grains in the interstellar medium. We make the first mea- 
surement of the neon ionization fraction in the ISM. We find Ne II/Ne Iw0.3 and Ne III/Ne I«0.07. 
These values are larger than is expected given the measured ionization of interstellar helium. For 
Ne IX, our results confirm the detection of the hot ionized interstellar medium of the Galaxy. 
Subject headings: ISM: general — X-rays: ISM — X-rays: binaries 



1. INTRODUCTION 

High-resolution X-ray spectral observations from the 
Chandra X-ray Observatory and XMM-Newton offer a 
new tool with which to study the interstellar medium 
(ISM). The ISM affects X-ray spectra in two ways: pho- 
toelectric absorption, particularly at low energies (0.1- 
10 keV), and scattering by dust grains, producing X-ray 
halos. At X-ray energies, absorption features are primar- 
ily from the excitation and ionization of inner-shell (K- 
shell) electrons, although for high-Z elements like iron, 
L-shell absorption edges are also detectable. Measure- 
ments of the ISM absorption features in the spectra of 
X-ray binaries allow us to study the abundances and ion- 
ization fractions for a large number of elements, analo- 
gous to optical and ultraviolet observations of stars that 
measure ISM absorption features. Additionally, high- 
resolution X-ray spectral measurements provide a new 
method for determining the elemental depletion of the 
ISM and composition of interstellar dust. 

The first attempt to measure ISM absorption edges in 
the X-ray s used the Einstein Focal Plane Crystal Spec- 
trometer ( Schattcn burg fc Canizaresll986j) and found ev- 
idence for the O I ls-2p (Ka) absorption resonance 
and a possible O II edge. Recent Chandra and XMM 
results have revealed more comple x structure around 
the oxygen J f-she ll absorption edge (|Paerels_et ii alj|20(^ 
[ Schulz et al.ll2002HTakei et al.ll2002Hde Vries etld1l200.4 
Uuett et alJ I2004L hereafter Paper I). All of the high- 
resolution work points to a complex system of absorp- 
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tion features at or near the oxygen edge, which is not the 
simple step function used in absorption models. Similar 
high-resolution structure is expected at all of the absorp- 
tion edges. 

Previous results from Einstein and EXOSAT grat- 
ing observations of the bright low mass X-ray binary 
(LMXB) Sco X-l point to an underabundance of oxygen 
and an overabund ance of nitrogen along the line of sigh t 
to the source (Kah n et al.lll984t iBrinkman et al.lll985tl 
Simila rly, using the abun dances of [A nders fc Grevessa 
l(T989l O/H = 8.51 x 1Q- 4 V IWeisskonf et a,l.l<|2004H found 
that the spectrum of the Crab pulsar was best-fit when 
oxygen was under abundant. When the newer abun- 
dances of IWilms et aTl IpiOfl O/H = 4.90 x 10~ 4 ) were 
used, oxygen was not required to be underabundant, al- 
though a slightly lower than nominal oxygen abundance 
was preferred, O/H = (3.33 ± 0.25) x 10~ 4 . We note 
that the measurement of the oxygen abundance relative 
to hydrogen is misleading since the X-ray band can not 
directly measure the hydrogen column density. Instead, 
the measurement is relative to the absorption from other 
metals which are converted to an equivalen t hydrogen 
colum n density using standard abundances. ITakei et alJ 
( 2002) compared the oxygen and neon column densities 
measured in the X-ray to the hydrogen column density 
from radio observations along the line of sight to Cyg 
X-2. Both the neon and oxygen abu ndances were found 
to be ~ 3/4 the solar abundances of lAnders fc Grevessd 
( 1989) but i n line with newer e stimates of ISM abun- 
dance s (e.g., IWilms et all l2000() . Finally, lUeda et all 
(2005) found an overabundance relative to oxygen for 
magnesium, silicon, sulfur, and possibly iron (a ssuming 
solar abundance ratios from I Anders fc Grevess e 1989) in 
the ISM absorption along the line of sight to three bright 
neutron star LMXBs. While many X-ray binaries have 
absorption which is well described by standard abun- 
dances, some systems have shown unusual abundance ra- 
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Fig. 1. — Left panel: Iron L-shell cross-section as a function of wavelength for metallic iron IKortright fc Kimll2000l) . The positions of 
the L,2 and L3 edges are marked. Right panel: Neon if-shell cross-section as a function of wavelength JGorczvca 2000). The positions of 
the ls-3p, ls-4p, and ls-5p transitions are marked. Only the ls-3p transition is resolved by Chandra. 



tios whic h have been attribu t ed to material local to these 
systems ilPaerels et al.ll20011 ilJuett et al.ll200lllLee et alJ 
120021 Uuett fc Chakrabartvl 120031) . Therefore, an accu- 
rate measurement of the relative abundances along many 
lines of sight is necessary in order to interpret data for 
sources which show unusual abundance ratios. 

We also mention a new result that found a significantly 
larger Ne/O ratio than in standard solar or ISM mod- 
els. U sing the X-ray emissio n lines from nearby solar-like 
stars, IDmke Test^l ll200^ found an average Ne/O ra- 
tio of 0.41, or 2.7 times greater than the standard ISM 
or solar abundances. If the Sun has a similar Ne/O ratio 
then this would solve th e solar helioseismology problem 
fsee iBahcall et all 120051) . We note that this work has 
sparked some controversy and claim s of smaller Ne/O ra- 
tios in the Sun were quickly made ijSchmelz et alJl2005l 
lAsnlund et all l2005aj) . In either case, it is not readily 
apparent that solar coronal measurements should equal 
the ISM abundances. 

It was recognized early in the history of X-ray spec- 
troscopy that the composition of dust grains could be 
probe d by spectros copic measurements of absorption 
edges l)Martinlll970(K Absorption edges from molecules 
include X-ray absorption fine structure (XAFS) due to 
the influenc e of nearby atoms on the absorption process 
ijWool Il995j) . Such features would distinguish between 
the molecular and atomic contributions to the edge struc- 
ture. In addition, the resonance structure of the absorp- 
tion edges is expected to be shifted, relative to the atomic 
position, in molecular absorption edges due to the molec- 
ular binding energy. The edges from magnesium, silicon, 
and iron are expected to be dominated by contributions 
from molecular, as opposed to atomic, forms given the 
high levels of depletion measured (Xd us t/^totai ^ 95%), 
while carbon, oxygen, and ni trogen are only mildly de - 
pleted ( A"d U st / A'totai ^ 60%: iSavaee fc Sembachlll996ft . 
Several authors have cl aimed detection of the ls-2p tran- 
sition in iron oxi des ilPaerels et alJ 120011: iSchulz et all 
120021 iTakei et al.i r20021. although we showed that this 
resonance coincides with the ls-2p transition from O II 



makin g the identificati on of iron oxides debatable (Pa- 
per iL lLee et alJ ()2002j) claimed a tentative detection of 
XAFS in the silicon absorpt ion edge of the X-r ay binary 
GRS 1915+105. Similarly II led a et all p)05j) claimed 
that the silicon and magnesium edges in the spectra of 
three LMXBs show XAFS consistent with the expected 
structure from silicates. In addition, depletion of ele- 
ments in grains can reduce the effect ive cross-section fo r 
absorption due to screening (see e.g.. lWilms eTalE oOO). 
This effect will lead to deviations in the measured abun- 
dances from the true interstellar values. 

Along with neutral and molecular absorption edges, 
features from ionized species are also expected. We 
showed that low ionization lines from oxygen are found 
in the spectra of a number of X-ray binaries (Paper I). 
From these lines, a measurement of the large -scale ioniza- 
tion fra ction of oxygen in the ISM is made. iMiller et all 
( 2004a) identified Ne II and Ne III absorption lines in the 
spectra of two black-hole X-ray binaries. They attributed 
these to material local to the systems, although such lines 
may arise in the ISM similar to the oxygen lines. Highly 
ionized lines may also be expected from t he hot ionized 
phase of the ISM. iFutamoto et all l)2004|) claim the de- 
tection of an O VII absorption line in the spectrum of 
the LMXB 4U 1820-30, which they attribute to the hot 
ISM. Recently, lYao fc Wanel l)2005|) detected Ne IX ab- 
sorption in the spectra of a number of LMXBs. Low- 
ionization lines can be used to better understand the 
relationship between the warm ionized and warm neu- 
tral phases of the ISM. The high-ionization features in 
Galactic LMXBs probe the hot ionized phase of the ISM 
and have direct implications for the detection of z = 
highly ionized oxygen and neon absorption features in 
the spectra of quasars, which have been claimed to be 
due t o the local group w arm-hot intergalactic medium 
(e.g.. iNicastro et al.ll2003 |) or al ternatively a hot Galac- 
tic corona (JSembach ct al. 2003(). 

In this paper, we study the absorption edges of neon 
(if-shell) and iron (L-shell) in nine bright X-ray binaries. 
These edges are most prominent in systems with Galactic 
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hydrogen column densities of 10 21 -10 22 cm -2 . Many of 
the X-ray binaries used in this work were also part of the 
oxygen edge study, allowing us to measure the relative 
abundances of oxygen, iron, and neon. These elements 
also sample a variety of expected depletion levels: neon 
should have no depletion, oxygen is mildly depleted, and 
iron is heavily depleted. We will look for signatures of 
depletion by comparing the abundance ratios of the ele- 
ments. Finally, we will measure the ionization of neon, 
using the identified low- and high-ionization absorption 
features. 

2. EXPECTED STRUCTURE OF THE IRON AND NEON 
ABSORPTION EDGES 

In Paper I, we presented an overview of the basic 
atomic physics necessary to understand LT-shell absorp- 
tion in the X-ray band. We defer the reader to the dis- 
cussion presented there and focus in the following specif- 
ically on the iron L-shell and neon i^-shell absorption 
edges. 

The iron edge located at 17.5 A is due to the exci- 
tation and ionization of L-shell (n = 2) electrons. The 
high-resolution structure of the edge is resolved into three 
distinct edges, L\, L 2 , and L 3 . The L\ edge at 14.7 A, 
due to transitions from the 2s level, is the weakest of the 
edges and we neglect it in our study since our data are 
not sensitive enough to measure this absorption. The L 2 
and L3 edges, due to transitions from the 2pi/ 2 and 2p 3 / 2 
levels, are located at 17.2 and 17.5 A, respectively (see 
Table [T] and Figure CJ. 

To model the cross-sectio n from Fe-L, we use the m ea- 
surement of metallic iron bvlKortright fe Kiml lj2000j) . A 
comparison of the lKortright fc Kiml l|2000) cross-section 
with lower-resolution determinatio ns ijHenke et alJll993l 
iVerner et alJll993t lChantlerlll995|) finds that the values 
compare to within 15% at a wavelength of 16.8 A, where 
the high-resolution structure is no longer a factor . Other 
analyses have found that the iKortright fc K im ( 2000) 
measurement reproduces the high-res olution structure 
of the iron edg e quite well (e.g., iLee et all 120011 
ISchulz et al.ll2002|K 

We expect the cross-section from iron in dust to differ 
in two important aspects. First, the shape and position 
of the edge can differ depending on the binding of iron 
within a molecule. The two most common forms of iron 
in molecules are ferrous (Fe 2+ ) and ferric (Fe 3+ ). Spec- 
troscopic studies of the Fe-L edge reveal that the posi- 
tion of the L3 edge varies from 17.517 A for ferrous iro n 
to 17.475 A for ferric iron <)van Aken fc Liebscherl l2002'). 
Molecules with a mix of the two forms show a double- 
peaked profile. The L? edge positio n and structure also 
varies (jvan Aken fc Liebscherll2002[). Olivines, a favored 
form of interstellar dust (see e.g. . lDrainell2003|) . contains 
pure ferrous iron. Iron oxides can contain ferrous (FeO), 
ferric (Fe20s), or ferrous ferric iron (Fe304), depending 
on the composition. Therefore, determination of the po- 
sition and shape of the L3 edge can yield information on 
the composition of interstellar dust. 

The second difference is the strength of the cross- 
section for absorption. In grains, the optical depth can 
be large (~1), therefore absorption will occur primarily 
on the surface. Atoms located in the inner regions of 
the grains are shielded by the outer layers. This pro- 



Table 1. Iron and Neon Edge and Line Positions 



Predicted Measurcd a 
Feature A (A) A (A) 

Kortright fc Kim 2000 

Fe-L 3 edge b 17.537 17.498±0.003 

Fe-L 2 edge b 17 - 2 . 26 17.188±0.003 

van Aken & Licbscher 2002 
Fe^+-L 3 edge b 17.517 17.498±0.003 
Fe 2 +-L 2 edge b 17.206 17.188±0.003 
Fe 3 +-L 3 edgc b 17.475 17.498±0.003 
Fe 3 +-L 2 edgc b 17.156 17.188±0.003 

Gorczyca 2000 

Ne 1 ls-3p c 14.298 14.295±0.003 

Nc I ls-ip 14.273 

Ne I ls-5p 14.263 

~ Gorczyca fc McLaughlin 2005 

Ne 11 ls-2p 14.605 14.608±0.002 

NeIIIls-2p 14.518 14.508±0.002 

Behar fc Netzer 2002 
Nc 11 ls-2p 14.631 14.608±0.002 

NeIIIls-2p 14.526 14.508±0.002 

Ne IX ls-2p 13.448 13.4439±0.0013 

a Measured in this work. 

b Position measured at wavelength of maximum absorption. 
c Position of ls-3p transition will be coincident with edge wave- 
length when using a standard edge model. 



duces a net reduction in the cross-section for elements 
located in dust grains. We have estimated this effec t 
for iron, following the description in Wil ms et alJ {2000), 
and find that the cross-section for absorption of iron in 
dust grains is 88% of the cross-section for gaseous iron. 
We note that this calculation requires a number of as- 
sumptions, including the relative abundances and deple- 
tions of the interstellar elements (see Appendix [SJ . A 
more detailed study of these assumptions is presented in 
Section IP 

The neon if-shell (n = 1) absorption edge is located 
at 14.3 A ijGorczvcal I2000K The neon edge, like oxy- 
gen, is composed of a series of resonance lines end- 
ing at the series limit or ionization energy. Neon, and 
higher Z elements, do not have a ls-2p transition due 
to a closed 2p s hell in the neutr al atomic ground state 
configuration. IGorczvcal lj2000[) calculated the cross- 
section for neutral neon absorption and found that his 
results closely matched experimental determinations of 
the cross-section (see Figure QJ. The resolution of Chan- 
dra and XMM is near the limit needed to resolve the 
ls-3p resonance line from the rest of the edge structure 
of neon, therefore we expect that when using standard 
edge models the measured edge positions will be consis- 
tent with the position of the ls-3p resonance (see Ta- 
bled). 

We also expect ionized neon to be present in the 
ISM. The primary transition we will find is from 
the ls-2p transition in ionized neon. Unlike neutral 
neon, the \s-2p t ransit ion is available for ionized forms. 
iBehar fc Netzerl l)2002l) calculated the wavelengths and 
oscillator strengths of the strongest ls-np (primarily ls- 
2p) transitions for a variety of elements and ionizations. 
The authors estimate that the errors on the calculated 
values are <0.2% for wavelengths and <10% for oscilla- 
tor strengths and autoionization rates. They also note 



4 



JUETT ET AL. 

Table 2. Observation Log 



Source Name ObsID Observation Date Time (ks) I b Dist. (kpc) \z\ (pc) Nq (10 18 cm~ 2 ) a 



4U 1820—30 


1021 


2001-07-21 


9 


2.79 


—7.91 


7.6 b 


1050 


1 -° -0.14 




1022 


2001-09-12 


9 










GX 9+9 


703 


2000-08-22 


21 


8.51 


+9.04 






q r>+1.2 

2 -0.6 


Ser X-l 


700 


2000-06-13 


76 


36.12 


+4.84 


9.5-12. 7 C 


800-1070 


Cyg X-l 


107 
3407 
3724 


1999-10-19 

2001- 10-28 

2002- 07-30 


9 
21 
14 


71.33 


+3.07 


2.0 d 


110 


6.3+1.4 
2.4+0.3 
2.9+0.3 


Cyg X-2 


1102 


1999-09-23 


29 


87.33 


-11.32 


11. 4-15. 3 C 


2240-3000 


1.1+0.2 


4U 1636-53 


105 


1999-10-20 


29 


332.91 


-4.82 


3.7-4.9 c 


310-410 


2 6+ 10 




1939 


2001-03-28 


26 










GX 339-4 


4420 


2003-03-17 


74 


338.94 


-4.33 


>6 C 


>450 




4U 1735-44 


704 


2000-06-09 


24 


346.05 


-6.99 


8.0-10.8 c 


970-1310 


3.4+1.2 


GX 349+2 


715 
3354 


2000-03-27 
2002-04-09 


9 
26 


349.10 


+2.75 









a Oxygen column density measurement from lJuett et all (2004). 
'iKuulker s et al.| 12003. and references therein). 
c Jonkcr & Nclcmans 12004, and references therein). 
IGierlinski et al.l 119991) . 



that the accuracy increases for higher ionization states. 
We have also calculated photoelectric absorption cross- 
sections for Ne I I and Ne HI using the same m ethod as 
for neutral neon ijGorczvca fc McLaughlinl20 05) . We list 
the predicted positions of the Ne II, Ne III, and Ne IX 
lines from both sources in Tabled Comparing the pre- 
dicted wavelengths for the Ne II \s-2p transition, we find 
that the two calculations differ by ~25 mA, on the or- 
der of the error given bv iBehar fc Netzerl (J2002) . The 
oscillator strengths are consistent within a few percent. 
We note that the pub lished natural width for Ne II in 
IBehar fc Netzerl J2002) is incorrect and that the correct 
value is 2.710xl0 14 s -1 (Behar 2005, priv. comm.). 

Using the full cross-sections for the neon species, we 
have checked that the neutral neon edge does not over- 
lap with transitions from Ne II and Ne III. The Ne II 
resonance series limit is at «13.7 A while the Ne III limit 
is at «13.3 A. No resonance lines from Ne II or Ne III 
are coincident with the neutral neon edge. The edges 
are easily distinguishable given the spectral resolution of 
Chandra. 

3. OBSERVATIONS AND SPECTRAL FITTING 
PROCEDURE 

Our study was based on archival observations of 
nine bright X-ray binaries obtained with Chandra us- 
ing the High Energy Transmission Grating Spectrometer 
(HETGS) in combination with the Advanced CCD Imag- 
ing Spectrometer fACIS: iCanizares et all 12005ft . The 
HETGS spectra are imaged by ACIS, an array of six 
CCD detectors. The HETGS/ACIS combination pro- 
vides both an undispersed (zeroth order) image and dis- 
persed spectra from the gratings. The various orders 
are sorted using the intrinsic energy resolution of the 
ACIS CCDs. The HETGS carries two types of trans- 
mission gratings: the Medium Energy Gratings (MEGs) 
with a range of 2.5-31 A (0.4-5.0 keV) and the High 
Energy Gratings (HEGs) with a range of 1.2-15 A (0.8- 
10.0 keV). The first-order MEG (HEG) spectrum has a 
spectral resolution of AA = 0.023 A (0.012 A). The abso- 
lute wavelength accuracy of the MEG (HEG) is ±0.011 A 



(0.006 A). 

In Table [3 we list the observations used in this analy- 
sis. The data sets were reduced using the standard CIAO 
threads 5 . In some cases, the zeroth-order data was not 
telemetered. For these observations, we estimated the 
zeroth-order position by finding the intersection of the 
grating arms and readout streak. After spectral extrac- 
tion, the accuracy of the estimated zeroth-order position 
was verified by comparing the wavelengths of strong in- 
strumental edges in both plus and minus sides of the 
spectra. Events collected during thermonuclear X-ray 
bursts were excluded from this analysis. Since the in- 
dividual observations of 4U 1820-30 and 4U 1636-53 
were brief, the spectra from multiple observations were 
combined for each source in order to improve statistics. 
Detector response files (ARFs and RMFs) were created 
for the +1 and -1 MEG and HEG spectra. The ARFs 
include the standard correction for the time-dependent 
change in the response due to a contaminant on ACIS. 
We then added the +1 and —1 orders to produce a com- 
bined MEG and HEG spectra from each source. Back- 
ground spectra were also extracted. The MEG data were 
used for both the iron and neon edge fits, while the HEG 
data were used for the neon edge when reasonable signal- 
to-noise was found. All fitting was performed in count 
space. In this paper, we present the flux-corrected spec- 
tra in order that the reader may better visualize the inter- 
stellar features. These flux spectra arc created in ISIS us- 
ing the tool f lux_corr which yields model-independent 
flux spectra. 

3.1. Iron Edge Fit 

We fit only the 16-18.5 A wavelength range and 
binned the data to ensure at least 10 counts per bin. 
An absorbed power-law continuum model was used for 
all sources. We allowed the power-law model parame- 
ters to vary during the fi ts. The tbvarabs absorption 
model (|Wilms et aLl EoOOl was used with the iron abun- 
dance set to zero to allow for an explicit modeling of 
the iron absorption. The equivalent hydrogen column 

5 http:/ /asc. harvard.edu/ciao/threads/ 
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Table 3. Fe Edge Parameters 



Source 


A (A) 


N Fc (10 16 cm" 2 ) 


/itt i son 3n 


1 ( .oUOztU.Ul t 


z< 1 4-1 O 


GX 9+9 


17.512±0.015 


7±2 


Ser X-l 


17.502±0.010 


15±3 


Cyg X-l ObsID 107 


17.495±0.009 


14.7 +2 | 


Cyg X-l ObsID 3407 


17.501±0.007 


lo.etfl 


Cyg X-l ObsID 3724 


17.491±0.004 


9.6±0.9 


Cyg X-2 


17.510±0.012 


5.6±1.4 


4U 1636-53 


17.495±0.010 


10±2 


GX 339-4 


17.505±0.007 


14 4+ 2 ' 2 


4U 1735-44 


17.51±0.02 


6±3 


GX 349+2 ObsID 715 


17.52±0.04 




GX 349+2 ObsID 3354 


17.49±0.02 


26lf 



Table 4. Ne Edge Parameters 



Source 


A (A) 


\r /-i rtl7 — 2 \ 

N Ne (10 1 ' cm z ) 


4U 1820-30 


14.294±0.012 


3.3±1.2 


GX 9+9 


14 304+ - 008 
J-*-3U*_0 .013 


3.9±1.3 


Ser X-l 


14.293±0.007 


7 67+ - 19 

'■ ol -0.7 


Cyg X-l ObsID 107 


14.317±0.007 


7.1±1.0 


Cyg X-l ObsID 3407 


14.293±0.004 


'• 4 -0.3 


Cyg X-l ObsID 3724 


14.256±0.012 


8.6±0.7 


Cyg X-2 


14.294±0.010 


9 O+0.9 
i *°— 0.3 


4U 1636-53 


14.291±0.012 


4.9±1.1 


GX 339-4 


14.281±0.012 


6.08+.0.13 


4U 1735-44 


14.315±0.013 


n 9 + 1.3 
z ' z -0.8 


GX 349+2 ObsID 715 


14.30±° ;gi 




GX 349+2 ObsID 3354 


14.284l°;gf 3 


13±3 



density (-ZVh) of the tbvarabs model was set to equal 
a multiplicative factor times the iron column density 
(iVp e ), where the multiplicati ve factor was the in verse 
of the iron ISM abundance of iWilms e~ al. ( 2000). We 
modeled the iron absorption edge by using a custom 
multiplicative model ba sed on the iron cross-section of 
IKortright fc Kiml l|2000D . The model parameters were 
Np c and the edge wavelength, defined as the position 
of the maximum cross-section of L3 edge. Both were 
allowed to vary in the fits. The best-fit N-p e and edge 
wavelengths are given in Table |31 All error estimates are 
90%-confidence levels, unless otherwise noted, and do not 
include instrumental errors. Figures El and |21 show the 
data and best-fit model for each spectrum. 

3.2. Neon Edge Fit 

We fit only the 13.1-15.4 A wavelength range and again 
binned the data to ensure at least 10 counts per bin. The 
continuum modeling used the same procedure as for iron 
(see §|OJ. We modeled the neutral neon absorption edge 
by using a cust om multiplicative model based on the neon 
cross-section of lGorczvcal l)200fl) . The model parameters 
were JVnc and the edge wavelength, defined as the posi- 
tion of the ls-3p resonance. Both were allowed to vary 
in the fits. Also included were three Gaussian absorp- 
tion lines to measure the absorption by Ne II, Ne III, and 
Ne IX. Initial fits fixe d the position of the se lines to those 
found for GX 339-4 ijMiller et alJl2004aj) and the widths 
(sigma) to 0.005 A (the instrument resolution). The ear- 
lier GX 339—4 result represents the best position identifi- 
cation for the low-ionization lines and was therefore used 
as the reference. Our results confirm the position of these 
lines. For sources where all three absorption lines were 
significantly detected, as opposed to only upper limits on 
their flux, we then refit the data, allowing the line posi- 
tions and widths to vary. The best-fit parameters were 
used to determine the equivalent widths (EWs) of the 
lines. The best-fit N-^ e and edge wavelengths are given 
in Table 01 while the best-fit line positions, widths, and 
EWs are given in Tabled Figures 01 and show the data 
and best-fit model for each spectrum. 

4. RESULTS AND ANALYSIS 

We now compare the results of our analysis for the 
sample of sources. From this we will determine the line 
and edge positions, the relative abundances of the neutral 



Table 5. Ne line parameters 


Source 


A (A) 


a (A) 


EW (mA) 


4U 1820-30 


14.5065 (fixed) 


0.005 (fixed) 


> -5 




14.6067 (fixed) 


0.005 (fixed) 


> -7 




13.4420 (fixed) 


0.005 (fixed) 


-4 8 +1 n 9 n 

^'°-2.9 


GX 9+9 


14.511±0.013 


0.013±0.010 


" 9 -7 




14.610±0.010 


<0.009 


-5±3 




13.438±0.008 


<0.010 


^+2 
-6I3 


Ser X-l 


14.510+.0.004 


<0.014 


-7-2-1.7 




14.606±0.006 


0.012±0.006 


-12±3 




13.446±0.002 


<0.007 


-7.7±1.4 


Cyg X-l 


14.5065 (fixed) 


0.005 (fixed) 


> -6 


ObsID 107 


14.6067 (fixed) 


0.005 (fixed) 


-7±3 




13.4420 (fixed) 


0.005 (fixed) 


> -3.6 


Cyg X-l 


14.511±0.005 


<0.013 


-4±2 


ObsID 3407 


14.612±0.004 


0.009±0.006 


-6±3 




13.4420 (fixed) 


0.005 (fixed) 


> -2 


Cyg X-l 


14.5065 (fixed) 


0.005 (fixed) 


-3. Oil. 9 


ObsID 3724 


14.6067 (fixed) 


0.005 (fixed) 


-3.2±1.9 




13.4420 (fixed) 


0.005 (fixed) 


-5.3±1.8 


Cyg X-2 


14.5065 (fixed) 


0.005 (fixed) 


> - 3.5 




14.6067 (fixed) 


0.005 (fixed) 


-6.2±2.7 




13.4420 (fixed) 


0.005 (fixed) 


> -5 


4U 1636-53 


14.502±0.008 


0.013±0.010 


-7±4 




14.601±0.012 


<0.04 


-10±7 




13.442±0.006 


<0.015 


-5±2 


GX 339-4 


14.506±0.003 


0.008±0.003 


-9.4±1.3 




14.607±0.003 


0.011±0.003 


-12±2 




13.4425±0.0019 


0.005±0.002 


-11.6±1.3 


4U 1735-44 


14.5065 (fixed) 


0.005 (fixed) 


> -6 




14.6067 (fixed) 


0.005 (fixed) 


> -12 




13.4420 (fixed) 


0.005 (fixed) 


'-4 


GX 349+2 


14.5065 (fixed) 


0.005 (fixed) 


> -13 


ObsID 715 


14.6067 (fixed) 


0.005 (fixed) 


> -37 




13.4420 (fixed) 


0.005 (fixed) 


> -20 


GX 349+2 


14.5065 (fixed) 


0.005 (fixed) 


> -21 


ObsID 3354 


14.6067 (fixed) 


0.005 (fixed) 


> -27 




13.4420 (fixed) 


0.005 (fixed) 


> -18 



and ionized species, and depletion fractions. 

Our data are well fit by t he metallic iron cross-section 
of IKortright fe Kiml (J2000), although requiring a small 
shift in the position of the L3 edge. This is compara- 
ble to previous studies o f the iron edge fe.g.. ILee et 
l200lHSchulz et al J 12002^ . The weighted mean position 
of the Fe-L3 edge (taken as the position of maximum 
optical depth) is 17.498±0.003 A (errors are 90% confi- 
dence and do not account for the absolute instrumental 
wavelength accuracy of ±0.011 A; see Figure |SJ) . Even 
including the absolute error, this value is significantly 
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Fig. 2. — Flux spectra of the iron L-shell absorption region for the seven neutron star LMXBs included in this study. The data have 
been binned by a factor of two for illustrative purposes. The dashed lines indicate the positions of the L2 and L3 edges. 



different from the measured positions of both the fer- 
rous and ferric edges and li es between the two values 
ijvan Aken fc Liebscherl 12002]) . This may suggest that 
the ISM contains a combination of ferrous and ferric 
forms of iron. But we caution that the energy resolu- 
tion of the laboratory data (0.8 eV=20 mA at 17.5 A; 
Ivan Aken fc Liebscherl l2002'l is large enough to make ei- 
ther value consistent with our measurement. The spacing 
of the Li and L3 edges is well fit by our model and is 
more consistent with ferrous iron. This would imply that 



the laboratory measurement of the wavelength of ferrous 
iron is off by ~19 mA. 

If both ferrous and ferric forms are present in the ISM, 
we would expect a double-peaked structure in the com- 
bined edge. We searched the residuals from the best-fit 
for our data. A negative residual at 17.40 A was found, 
although with low significance, < 3<r. Even if true, the 
wavelength of the feature is too small to be attributable 
to a contribution from ferric iron, which we would expect 
at 17.45 A if our best-fit position is from ferrous iron. In 
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Fig. 3. — Flux spectra of the iron L-shell absorption region for 
the two black hole X-ray binaries included in this study. The data 
have been binned by a factor of two for illustrative purposes. The 
dashed lines indicate the positions of the L2 and L3 edges. 



addition, a check of the effective area versus wavelength 
for the sources with the most prominent residuals found 
sharp variations in the effective area, a signature of hot 
pixels or columns, in the iron wavelength range. We also 
note the positive residuals just high of the edge position 
(17.6 A). We know of no atomic or molecular features 
that would produce such a residual. It is also possible 
that these residuals are due to hot pixels or columns. 

For the neutral neon edge, we find a weighted mean 
position of the ls-3p transition of 14.295±0.003 A (see 
Figure 01. This is consistent w ith laboratory measure- 
ments (e.g.. iCoreno et aflll99 9). We note that the the- 



oretical position of iGorczvcal l)2000j) was calibrated to 
this laboratory measurement. For the ionized neon lines, 
we find a weig hted mean position of 14.608±0.002 A for 
Ne II, 14.508±0.002 A for Ne III, and 13.4439±0.0013 A 
for Ne IX. The Ne IX position is consistent with theoret- 
ical and observational measurements, accounting for the 
absolute instrument error. The Ne II and Ne III lines, 
however, are «20 mA off of the theoretical predictions. 
This is within the quoted error s for the theoretical calcu- 
lations ijBehar fc Netzerll2002j) . In this case, the error is 
more significant for the theory than the observations. We 
encountered a similar effect when identifying the neutral 
and low-ionization lines from oxygen (Paper I) . We note 
that our best-fit positions for th e Ne II and Ne III li nes 
are consistent with the results of iMiller et al.l l)2004aj) . 

4.1. Relative Abundances of Oxygen, Iron, and Neon 

Using the column density measurements, including the 
oxygen column density measurements of Paper I, we can 
now directly measure the abundance ratios in the X-ray 
domain. This measurement is a function not only of the 
interstellar abundances, but also of the depletion of the 
elements into dust grains. We have chosen to use neon 
as the calibration point for these measurements, thereby 
providing O/Ne and Fe/Ne ratios, for two reasons. First, 
as a noble gas, neon should not be depleted. Second, as 
the lowest wavelength edge of the set presented here, it 
lies near the middle of the Chandra wavelength band. 
Observations of higher column density sources will show 
edges from neon, magnesium, and silicon. Therefore, 
abundance ratios relative to neon will allow us to di- 
rectly compare the ISM abundances of various elements, 
e.g., oxygen and silicon, which could not be done with a 
single source in the X-ray waveband. 

We calculated the abundance ratios (by number) for 
the entire sample and for the neutron star (NS) LMXBs 
only. While the black hole systems in our sample are 
known to show intrinsic features iSchulz et al.l I2002t 
iFeng et alJl20?)3 IMiller et al.ll2"f)03t IMiller et alJl20(CT 
the NS LMXBs show no emission or absorption features 
that are intrinsic to the system. NS LMXBs which do 
show intrinsic features are not included in this study. 
We therefore take the NS subsample to be indicative of 
the ISM properties and show the black hole systems rel- 
ative to this subsample. We note that the inclusion of 
the black hole systems in this study may seem unneces- 
sary given their intrinsic features, but they provide the 
highest signal-to-noise spectra allowing us to make more 
accurate measurements of the line and edge positions. 

In Figure [5] we compare the neon and oxygen column 
densities for our sample. For all the systems, we find a 
O/Ne ratio of 3.7±0.3. When we restrict our sample to 
the NS systems, we find 0/Ne=5.4±1.6. The difference 
in the measurements is driven by the low oxygen column 
measurements for two of the Cyg X-1 datasets. We note 
that while the oxygen column for Cyg X-1 shows a fac- 
tor of 2 difference between the three datasets, the neon 
column density measurements are all consistent within 
errors. This is likely due to material local to the source 
(see Paper I and references therein). 

Our abundance ratio measurement for the NS sys- 
tems is consist ent with 0/Ne=5.6 o btained from the ISM 
abundances of I Wilms et alJ l)2000|) . This is also consis- 
tent with standard solar abundance ratios (0/Ne=6.6; 
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Fig. 4. — Flux spectra of the neon ii"-shell absorption region for the seven neutron star LMXBs included in this study. The data have 
been binned by a factor of three for illustrative purposes. The dashed lines indicate the positions of the identified features. 



e.g.. lLodder3l2003UAsplund et alj |2005b). In terestingly 
our re sults exclude the abundance ratio of lDrake fc Testal 
( 2005), 0/Nc=2.4. One possible explanation is that our 
two studies cover a very different extent in the Galaxy. 
All of our sources are >2 kpc from the Sun, while 
the [Drake fc Testa ( 2005) sample of stars are all within 
100 pc. 

Our best-fit O/Ne ratio also allows for roughly 45% 
of the oxygen in the ISM to be located in grains (see 
derivation below and Appendix 0, although we note 



that the large error on the O/Ne ratio allows no real 
constraint. More precise measurements of the O/Ne ratio 
are needed to probe the depletion in the ISM. 

The Fe/Ne ratio, howe v er, sh ows a large devia- 
tion from the IWilms et alJ lj2000|) ISM abundance ra- 
tio of 0.309 (see Figure ITU|) . Including all systems, we 
find Fe/Ne=0.159±0.013, while for NS systems only, 
Fe/Ne=0.20±0.03. Both measurements are significantly 
different than the ISM value. The sense of the discrep- 
ancy is that we are finding less iron in the ISM than we 
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Fig. 5. — Flux spectra of the neon if -shell absorption region for 
the two black hole X-ray binaries included in this study. The data 
have been binned by a factor of three for illustrative purposes. The 
dashed lines indicate the positions of the identified features. 



would expect. The discrepancy between the measured 
and predicted Fe/Ne abundance ratio is likely the result 
of depletion in the ISM. As mentioned previously, atoms 
in dust grains have a lower effective cross-section than the 
same atoms in gas. Therefore, using the gas cross-section 
(as we have) would cause the column density measure- 
ments to be lower than the true value. 

Attributing the abundance ratio difference to deple- 
tion, we can estimate the amount of iron located in dust 
grains. The total optical depth of iron, tf c , is the sum of 
the contributions from gas and dust grains. The gaseous 
iron optical depth is the gas phase cross-section times 
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Fig. 6. — Best fit position of the iron L3 absorption edge, de- 
fined as the position of maximum cross-section. The weighted mean 
value of the line position is given. All quoted errors are statisti- 
cal only and do not include the instrumental wavelength error of 
±0.011 A (FWHM). 



the column density of iron in the gas phase, given by 
Cgas-^Fe,gas- The optical depth from iron in dust is the 
grain cross-section times the column density of iron in 
grains given by a gia ,i n N-pe, grain- The depletion factor, j3, 
measures the amount of the total iron column density 
found in grains, iVpe^rain = (3N-p c . Similarly, the gas col- 
umn density is N-p Bigas = (1 — (3)N-p c . Finally, the dust 
cross-section can be related to thegas cross-section by a 
factor, /, as shown in Appendix lAl Combining this we 
find 



5 (l-/3)^ Fc + /a g a S /3iVFc = <r gas iV Fe (l-/3+//?) 



Our analysis assumes that Tp c = <J sas N-p c . 



Therefore we 



can relate our measured iron column density, A^F e ,meas, to 



the true column density by iVpe.meas = ATp e (l — (3 + //?). 
Since we do not know the true N-p e , we must use abun- 
dances ratios to measure depletion. Using abundance 
relative to neon allows us to avoid additional depletion 
effects that would be present with oxygen. The relation- 
ship between the measured Fe/Ne abundance ratio and 
the ISM ratio is given by 



(Fe/Ne) r 



(1 -0 + f0) x (Fe/Ne) IS M , (2) 



assuming (3 = for neon. Given our measured Fe/N e 
ratio, the ISM abundance ratio from lWilms et alJ J2000), 
and the value of / calculated in Appendix^] we find that 
(3 = 2.9 ± 0.8, which is unphysical since (3 is constrained 
to the range 0-1. 

Obviously, one or more of our assumptions is invalid. 
It is possible tha t the ISM abundan ces of either iron or 
neon is incorrect. iWilms et all <f2000^ claim errors for the 
ISM abundances of 0.1 dex or higher. A neon abundance 
0.1 dex larger is enough to account for the discrepancy 
and constrain (3 > 0.5. A slightly greater neon abun- 
dance is not at present problematic for our measured 
O/Ne ratio, given the large error bars. Similarly, a lower 
iron abundance could also account for the difference. 

We have also tested the possibility that the dust 
grai n cross-section cou ld be misestimated. Assuming 
the IWilms et alJ (UpOQJ) abundances, for various deple- 
tion factors, /3=0.7-0.95, we find that Equation [5] gives 
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Fig. 7. — Best fit position of the ls-3p transition from Ne I. The 
weighted mean value of the line position is given. All quoted errors 
are statistical only and do not include the instrumental wavelength 
error of ±0.011 A (FWHM). 



not the average value over the edge which would give 

a lower cross-sectio n. In addition, the cross-section 

of iKortright fe Kiml ( 2000) compares well with lower- 
resolution measurements at wavelengths away from the 
high-resolution structure. Therefore the cross-section 
difference between the L3 maximum and the lower wave- 
length absorpt ion would have to be sub stantially greater 
than found bv IKortright fc Kiml l)2000|) . Since our data 
fit both the high-resolution structure and the lower wave- 
length absorption, we would expect to see deviations in 
our data from the chosen model if this was the case. Such 
deviations are not found and we therefore conclude that 
the iron edge model is accurate in its structure and scale 
to at least 10%. 

Finally, we address a final assumption in the calcula- 
tion of /, which i s a homogeneous distribution of the ele- 
ments in grains. iSpitzer fc Fitzpatrickl (|lQ93f) suggested 
that interstellar silicate grains may have a iron-rich core 
with a magnesium-mantle. Such a scenario would pro- 
vide more shielding to the iron atoms, producing a lower 
value of /. 
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Fig. 8. — Best fit positions of the ls-2p transitions from Ne II 
(top panel) and Ne III (bottom panel). All quoted errors are sta- 
tistical only and do not include the instrumental wavelength error 
of ±0.011 A (FWHM). 



values of /=0.4-0.7, significantly lower than our calcu- 
lated value. We have calculated /, as shown in Ap- 
pendix lAl for various assumptions, including the dus t 
model l)Mathis et al.lll977tlWeingartner fc Draindl2001fl . 
abundances, variations in depletion factors of the astro- 
physically abundant elements, grain mass density, and 
the gas cross-section for iron. Most of these parameters 
do not make a large difference to the value of /. 

Only the cross-section and grain mass density are able 
to make a significant effect, although they must be more 
than 5 times larger to bring the measured and pre- 
dicted results into agreement. A difference this large 
in the iron cross-section is unlikely. First, we are al- 
ready using the maximum value of the cross-section, 



The diffuse ISM has four major phases: cold neu- 

tral, warm neutra l, warm ionized, and hot ionized 

ijrleiles fc Kulkarni 11987ft . The strength of the ion- 
ized neon lines can be used to constrain the proper- 
ties of these phases. The Ne II and Ne III lines are 
most likely associated with the warm ionized medium 
(WIM), while the neutral edges will be primarily from 
the warm neutral medium (WNM). From an ultravi- 
olet observation of t h e Gal actic halo star HD 93521, 
ISpitzer fc Fitzpatrickl <|1993f) found th at the WIM co- 
incid e d with the WNM. Ot her studies ijrlowk fc Savage! 
Il999t iRevnolds et alJ Il995f) have found similar results. 
This allows us to make a direct comparison of the col- 
umn densities from the neutral and low ionization neon 
lines from which the ionization of the ISM can be mea- 
sured. 

The ionized abundances can be estimated from a curve 
of growth analysis (see Paper I) . We used the theoretical 
cross-section of Gorczyca (2005; in prep.) to calculate 
the curve of growth for Ne II and Ne III assuming ther- 
mal velocities of 20, 100, and 200 km s -1 . Using the os- 
cillato r strengths and natural widths of iBehar fc Netzerl 
(2002) gives consistent results. In Figure ITT1 we plot the 
EW in the ionized lines versus the neutral neon column 
density for each source. Given an ionized abundance, 
A = Nx/N^ci, one can then overlay the curves of growth 
calculated from the theoretical line properties. The best- 
fit abundance is determined by comparing the theoret- 
ical curves at different values of A to the data points. 
We note that there is a degeneracy between the thermal 
velocity and the ionized abundance. Given the low sig- 
nificance of our EW measurements, we only estimate the 
ionized abundances. We find that Ne II/Ne I«0.3 and 
Ne lII/Ne I«0.07. 

Our analysis attributes the ionized neon lines to inter- 
stellar material. In our oxygen edge study, we argued 
that an interstellar interpretation is correct, since the 
EWs of the lines correlated better with the neutral col- 
umn density than with the source luminosity (see Paper 
I). The same trend is suggested by the ionized neon lines, 
although at much less significance. Obviously, more ac- 
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Fig. 9. — Comparison of the neutral oxygen and neon column 
densities for our sample of sources. The best fit O/Ne ratio is 
plotted for all sources (solid line) and the neutron star LMXBs 
only (dotted line) . We have overplotted the ISM O/Ne ratio of 
IWilms e ^alj 1200(1 dash ed lin e) and the stellar coronal measure- 
ment of lDra kc & Testa (2005), dashed-dotted line). The best fit 
value for the NS LMXBs in is good agreeme nt with the ISM value . 
None of our values are consistent with the Drake & Testa (2005) 
O/Ne ratio. 



curate measurements are required to verify this trend. 

It is interesting to note that the EWs of the low- 
ionization neon lines from GX 339—4 give a result in 
line with the NS LMXB trend. This suggests that the 
low and neutral absorption along the line of sight to GX 
339—4 is consistent with an interstellar origin. On the 
other hand, the EWs for the low-ionization neon lines 
from Cyg X-l do not follow the same trend. In the oxy- 
gen study, Cyg X-l showed large variations in the column 
densities of the neutral and ionized lines (Paper I). We 
attributed those changes to circumstellar material in the 
binary. The variations in the neon features are less sub- 
stantial, which may imply something interesting about 
the source of the ionizing flux. 

The Ne IX absorption should originate in the hot phase 
of the ISM. With a higher temperature and lower den- 
sity, the hot phase will have a larger scale height than 
the warm phases which contain the neutral and low ion- 
ization lines. Therefore, a dir ect compariso n is no t ap- 
propriate. Instead, we follow lYao fc Wand" <g005) and 
compare the Ne IX column densities as a function of z, 
the distance from the Galactic plane (see Figure IT2*)l . We 
overplot their best-fit disk distribution. The Ne IX col- 
umn densities were calculated from a cu rve of growth 
assum ing a gas temperature of 2.4 x 10 6 K l)Yao k, Wand 
2005). Our results compare reasonably well with their 
best-fit model, which is not surprising given the large 
overlap of sources. Again, we see variations in the Cyg X- 
1 (open triangles) column densities. The data point that 
deviates the most from the model is from GX 339—4 
(filled triangle). This confirms that a sizable contribu- 
tion to the Ne IX line is due to a local contribution, 
which is reasonable given the detection of other hydro- 
genic and helium-lik e lines in the spectrum of this source 
(Miller e t~aTll2004aH . 
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Fig. 10. — Comparison of the neutral iron and neon column den- 
sities for our sample of sources. The best fit Fe/Ne ratio is plotted 
for all sources (solid line) and the neutron star LMXBs only (dot- 
ted lin e). We have overplotted the ISM Fe/Ne ratio of IWilms et all 
(2000, dashed line). Our best fit Fe/Ne ratio is significantly lower 
than the ISM value. We attribute this to depletion of iron into 
dust grains in the ISM. 



5. DISCUSSION 

We have presented a global study of the neon if -shell 
and iron L-shell absorption in the ISM. Combined with 
our previous study of oxygen absorption (Paper I), we 
have determined the neutral abundances ratios between 
oxygen, iron, and neon. We find 0/Ne=5.4±1.6, in line 
with current ISM abundances. Our Fe/Ne abundance of 
0.20±0.03 is lower than the expected ISM value of 0.309, 
suggesting that iron is depleted in the ISM. We note that 
unless depletion is taken into account in X-ray absorption 
models, a subsolar iron abundance may be found. This 
subsolar iron abundance can be due only to the ISM and 
does not a priori signify an astrophysically interesting 
result related to the source properties. 

We also detected absorption lines from Ne II, Ne III, 
and Ne IX in a number of the sources. We attribute 
these lines to the ISM. A curve of growth analysis found 
that the large-scale ionization of neon in the ISM is 
Ne II/Ne I«0.3 and Ne III/Ne I«0.07. The implied col- 
umn den sities in the Ne IX li nes are consistent with the 
model of lYao k, Wand l)2005|) for all sources except the 
black hole LMXB GX 339-4. 

5.1. Ionization of the Warm ISM 

The ISM of the Galaxy is typically modeled as a mul- 
tiphase medium (see e.g.. lMcKee fc Ostrikerlll977j) . The 
phases have different temperatures, densities, and ion- 
izations. The following discussion focuses on two of the 
phases, the warm neutral and warm ionized mediums 
(WNM and WIM, respectively). These two phases have 
similar temperatures (T ~ 10 4 K) but different hydro- 
gen ionization fractions. In the WNM, the ionized hy- 
drogen fraction H H/ Htntai ~ 0.1, whil e in the WIM, 
H II/H to tai > 0.8 (see Ho wk et al J 1200.1 and references 
therein). 

Our implied ionization fractions for Ne II and Ne III 
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Fig. 12. — Ne IX column density as a function of source dis- 
tance above the plane. Column densities were calculated from the 
EWs assuming a gas temperature of 2.4 X 10 6 K. For systems with 
no known distance, we assume a distance of 8.5 kpc, while for 
GX 339—4 we used the lower limit on the distance of 6 kpc. Slight 
shifts in \z\ were applied for illustrative purposes . Overplotted is 
the best fit disk distribution o flYao fe"W ang ( 2005). The diamonds 
mark the NS LMXB systems, the open triangle marks the black 
hole binary Cyg X-l, and the filled triangle marks the black hole 
LMXB GX 339-4. Notice that the data for GX 339-4 is the most 
discrepant from the model, confirming a local contribution to the 
Ne IX line. 



Fig. 11. — Curve of growth for the Ne II and Ne III absorp- 
tion features. Plotted are the measured values of the neutral neon 
column density, N^j e , and the EW in the ionized neon absorption 
lines. The diamonds mark the NS LMXB systems, the open tri- 
angle marks the black hole binary Cyg X-l, and the filled triangle 
marks the black hole LMXB GX 339-4. Overlayed are the the- 
oretical predictions for the relationship between the column den- 
sity and EW in the lines. The solid, dotted, and dashed lines 
indicate the predicted EWs for velocity dispersion of 20, 100, and 
200 km s — , respectively. Also included are predicted EWs for var- 
ious abundances (A = N^/N^el) relative to neutral neon. From 
these figures, we can estimate the relative abundances for Ne II 
and Ne III relative to Ne I. 



are greater than we found for oxygen, O II/O IwO.l and 
O III/O I«0.05 (Paper I). With higher first ionization 
potentials, we would expect the ionization of neon to 
track that of helium in the ISM. Measurements of the 
He I A5876/Ha line intensity ratio suggest a helium ion- 
ization fraction n(He II)/n(He to tai) ^ 0.27 in regions 
where the hydrogen is pr imarily ionized (i.e. in the WIM; 
IRevnolds fc Tufta lT9951. This has been taken to imply 
that the helium is primarily neutral in the WIM, since 
the alternative explanation, that helium is primarily fully 
ionized, appears inconsistent with the low ionization of 
other elements, particularly sulfur and nitrogen. 

Our measured neon ionization ratios include the con- 
tributions from all phases of the ISM. In the following, 
we attribute the neutral and low-ionization features to 
the WNM and WIM only. These phases should be the 
dominant contributors. We would like to estimate the 
Ne II/Netot ratio in the WIM to compare with the he- 
lium result. To do this, we assume that the Ne II and 
Ne III lines are produced in the WIM only, while the 



neutral neon absorption includes contributions from both 
the WNM and WIM. Given the higher ionization poten- 
tial of neon compared to hydrogen, we believe this is 
valid based on the hydrogen ionization in each phase. 
We must make two other assumptions: first, the Ne I, 
Ne II, and Ne III are the dominant contributors to the 
total interstellar neon budget, and second, that a sig- 
nificant fraction (> 75%) of the neutral neon is found 
in the WNM, as opposed to the WIM. Combining these 
assumptions, we find 

/NeII\ 

\ Ne tot / WIM 

(Ne II/Ne I) tot 

Ne Ivvm/Ne I tot + (Ne II/Ne I) tot + (Ne III/Ne I) tot ' 

(3) 

This work directly measures (Ne II/Ne I) to t and 
(Ne III/Ne I) tot , and we are assuming Ne Iwnvi/Ne Itot ^5 
0.25. 

Our results then yield Ne II/Ne to t > 0.4 in the WIM. 
This is different from the helium results, although we 
note that the exact values of the ionization fraction of 
Ne II and Ne III are dependent on the assumed ther- 
mal velocity of the ISM, which may make the differ- 
ence between the measured ionizations of helium and 
neon less significant. In addition, our measurement of 
the Ne III/Ne I ratio seems incompatible with standard 
OB star models for the ionization of the WIM (see e.g., 
iSembach et al.ll200(Tfl . Extra heating has been studied to 
explain the properties of optical emission and absorptio n 
lines (e.g.. iSlavin et alJl200fi lElwert fc Dettmarll2005|) . 
but these results concentrate on the elements of interest 
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in the optical. Therefore, neon ionization fractions are 
not given. One possible explanation for our neon ion- 
ization results is that the lines of sight probed by the 
X-ray studies contain a larger WIM fraction relative to 
the WNM. 

X-ray studies probe different lines of sight than tra- 
ditional optical and ultraviolet work. Many of our sight 
lines are towards the center of the Galaxy and the Galac- 
tic bulge. Optical and ultraviolet studies concentrate on 
the Perseus arm of the Galaxy which is roughly anti- 
Galactic center (120° < I < 150°). While X-ray allows 
us to probe larger distances and cover a greater extent of 
the Galaxy, this also makes it difficult to reconcile X-ray 
results with other ISM studies. Future work will use mul- 
tiwavelength studies of X-ray binaries to compare X-ray 
ISM results with optical and ultraviolet studies. 

5.2. The Warm Absorber m GX 339~4 

iMiller et all (|2004aj) identified the Ne II, Ne III, and 
Ne IX lines in the spectrum of GX 339—4 and attributed 
them to an AGN-like warm absorber. Our results show 
that the Ne IX line shows a clear excess column density 
over the expected interstellar component. The measured 
Ne IX column density for GX 339-4 is -3.9 x 10 17 cm" 2 . 
For Galactic latitude of GX 339—4 and a distance of 
6 kpc, the ISM model predicts a column density of 
7.8 x 10 15 cur 2 . While the exact distance to GX 339-4 
is unknown, the model column density as a function 
of Galactic scale height reaches a maximum value of 
2.5 x 10 16 cm~ 2 . From these values, we estimate that 
only «2-6% of the measured Ne IX column density to 
GX 339—4 is due to the hot ISM. This seems reasonable 
given the detection of other highly ionized lines in the 
spectrum of the source. 

The properties of the Ne II and Ne III lines how- 



ever, are completely consistent with the other LMXBs in 
our sample. The measured line widths are all w20 mA 
off fro m the theoretical calculations of iBehar fc Netzerl 
(2002), within the quoted errors for the calculations. 
Similar wavelength shifts were also required in our study 
of the oxygen edge (Paper I). The line widths are also 
consistent with the other sources in our sample, and are 
more likely indicative of the resolution of the instrument 
rather than an astrophysical origin. Finally, the im- 
plied ionized abundances of Ne II and Ne III from the 
GX 339-4 data are consistent with the NS LMXBs. If 
these lines had a local component, we should see devia- 
tions from these relationships. We note th at variations 
in the EW of these lines were claimed by IMiller et alJ 
( 2004b) but we caution that comparisons between Chan- 
(ira/HETG and XMM/RGS spectra are suspect for nar- 
row absorption features due to the factor of two differ- 
ence in the spectral resolution between the instruments. 
Narrow features easily detectable by Chandra can be 
substantially washed out by the RGS resolution. We 
therefore claim that the low-ionization neon features in 
GX 339—4 are completely consistent with an interstellar 
origin. 
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APPENDIX 

DERIVATION OF CROSS-SECTION FOR DUST GRAINS 

The derivation of the cross-section for dust grains presented in IWilms et al.1 lj2000|) includes the contribution of all 
elements. For our data, we need the iron dust cross-section only. Therefore, we present a rederivation of the calculation 
of the dust cross-section for a single eleme nt, Z. 

We begin by considering Equation A7 of Wil ms et alJ l|2000|) : 



dn 



[a 



grams = A H £g J ^ V geom [l - CX P (- (ff) ( A) )]da , (Al) 

where Ah is the hydrogen column density, £ g is the number of grains per hydrogen atom, dn gr (a) / da is the grain 
size distribution assuming spherical grains of size a, c gcom = ira 2 for spherical grains, (a) is the average photoelectric 
absorption cross-section of the grain material and (N) is the column density in a grain. 

For our purposes, we want r gra i ns of element Z. Ah can be directly related to the element column density Nz by 
Ah = Nz/Az, where Az is the ISM abundance of element Z. Since we are considering a single element, (a) is just 
the gas phase cross-section of that element, cx gas . 

To calculate (N) for a single element, we follow IWilms et all (HxxJ, and assume that the grains have a homoge- 
neous chemical composition and are spherical with size a. We consider only the contribution of element Z to (A), 
remembering that element Z makes up only part of the total mass of the grains. We find 

(AO = ^ , (A2) 

where pz is the mass density of element Z in interstellar grains and \xz is the atomic mass. The mass density pz is 
related to the total grain mass density p by 

pz = p x AzPzpiz J X! A zPzHz , (A3) 
z 

where (3z is the depletion fraction of element Z . The calculation of pz takes into account the contribution of all other 
elements to the total dust content of the ISM. 
For a single element, the number of grains per atom £ g is given by 

/ f dn„ r (a) 4 , , . , 

£ g = AzPzpz / / V* x p z x -vra 3 x da . (A4) 

Substituting the above relationships into Equation IA1I we find 

at a f°° dn„ T (a) 2 . Ap z a I f dn Rl {a) 4 , 

rgrains = N z (3 Z ^z J g v ' no 1 [1 - cxp(-cr gas )]da J ^ w p z -7ra 3 rfa . (A5) 

For our analysis, we are interested in the determining the value of /, the factor which relates the dust cross-section 
to the gas-phase cross-section, cr gra in = /c gas . 

The grain optical depth r gra i n = cr g rainAz, gr ain where Az !gra j n is the column density of element Z found in grains. 
The grain column density is given by the total column density times the depletion factor, Pz, i.e. A^ !grain = fizNz- 
Therefore, we can substitute for the grain optical depth, T gra j ns = f& glis PzNz- Solving for /, we find 

r f°° dn gY (a) 2 4p z er gas / f dn gI (a) 3 

/ = ~. / , - a 1 - exp( —^a)]da / / & x a 4 x da , (A6) 

4pzCT gas Jo « a «W / J da 



after some simplification. 
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We can now use the above equation to calcul ate / under various assumptions. Our initial assumptions use the 
ISM abundances an d depletion factors given bv I Wilms et al.1 l|2000l) . We take p = 1 g cm -3 and dn gT (a)/da from 
iMathis et all l)1977|) . We note that (Tg as is a function of energy. 

For iron, we assume a single value for cr g a S , which is the maximum value of the cross-section at the L3 edge. Using 
the maximum value of <7 gas gives the lowest value of /. We find / = 0.88 for iron. 

Similarly, for oxygen we can calculate the value of /. The dominant X-shell edge for oxygen occurs at the ls-np 
series limit for the S = 3/2 final state (see Paper I). Using the cross-section value at this limit to find the value of / 
yields / = 0.90. 



